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ABSTRACT 

We present Hubble Space Telescope observations in the far UV of the ultraluminous X-ray source 
in NGC 6946 associated with the optical nebula MF 16. Both a point-like source coincident with the 
X-ray source and the surrounding nebula are detected in the FUV. The point source has a flux of 
5 x 10 -16 ergs -1 cm~ 2 A -1 and the nebula has a flux of 1.6 x 1CU 15 ergs" 1 cm" 2 A" 1 , quoted at 1533 
A and assuming an extinction of Ay — 1.54. Thus, MF 16 appears to host the first directly detected 
ultraluminous UV source (ULUV). The flux of the point-like source is consistent with a blackbody 
with T w 30, 000 K, possibly from a massive companion star, but this spectrum does not create 
sufficient ionizing radiation to produce the nebular Hen flux and a second, hotter emission component 
would be required. A multicolor disk blackbody spectrum truncated with an outer disk temperature 
of ~ 16, 000 K provides an adequate fit to the FUV, B, V, I, and Hen fluxes and can produce the 
needed ionizing radiation. Additional observations are required to determine the physical nature of 
the source. 

Subject headings: black hole physics — galaxies: individual(NGC 6946) — ultraviolet:stars — X-rays: 
binaries 



1. INTRODUCTION 

Ultraluminous X-ray sources (ULXs) are bright, ir- 
regularly variable, non-nuclear, X-ray sources with 
apparent luminosities exceeding the Eddington limit 
for a 20Mq compact object. ULXs could represent 
a new class of black holes with masses intermedi- 
ate between stellar-ma s s and s upermassive black holes 
dColbert fc Mushotzkvl Il999t iMakishima et all 120001: 
iKaaret et all I2001D or may be stellar-mass black holes 
with unusually high accretion rates and/or beamed emis- 
sion. Over the past few years, optical and radio nebu- 
lae have been found which are s patially coinciden t and 
likely physically related to ULXs (iPakull fc Mirionill2002l; 
Roberts fc Colbertl [200H: IKaaret et all 12004 lLang et al l 



20071 ). Study of nebulae associated with ULXs may en- 



able us to determine the total radiation and particle 
fluxes from the ULXs if the nebulae are continuously 
powered by photoionization or particle flows. 

The nebula MF 16 in the nearby (5.1 Mpc) and 
nearly face-on spiral galaxy NGC 6946 has been de- 
scribe d as an 'ultraluminous super nova remnant com- 
plex' ijBlair. Fesen. fc Schlegel 120011 ) . The object is un- 
usually luminous for a SNR in optical line emission, e.g. 
the Ha l umin osi ty is 2 x 10 38 ergs ~ 1 , an d in X-rays 
(|Schlegell[l99l . IRoberts fc Colbert! (12001) found that 
the X-ray source is variable and, thus, not diffuse X- 
ray emission associated with a SNR, but rather emission 
from an X-ray binary containing a compact object. MF 
16 which had been previously thought of as a SNR is, 
instead, a nebula surrounding an X-ray binary. 

MF 16 shows many lines consistent with shock emis- 
sion. However, several of the lines (e.g. [O III], 
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Ha, and [N ii]) contain a narrow component, gener- 
ally associated with photoion ization instead of shocks. 
IBlair. Fesen. fc Sch legel (2001) concluded that unusually 
strong photoionization is needed, in addition to shocks, 
to explain the spectrum, but did not determine the 
source of ionizing photons. They also reported detec- 
tion of Hell line emission, but did not discuss the origin 
o f the line. 

lAbolmasov et all (|2008[) studied the optical emission 
from MF 16 with particular attention to the Hen A4686 
line. Their modeling shows that most of the power in 
the optical line emission comes from photoionization, 
although a component from shocks is still important. 
They find that an ultraviolet source more luminous than 
the observed ULX is needed to power the emission and 
that the UV source has a blackbody temperature of 
10 5.i5±o.io K and a i um i nos ity 1.2 x 10 40 ergs _1 . The 
luminosity is a lower bound on the source luminosity if 
the photoionized nebula is dens ity bounded. 

The UV source proposed by lAbolmasov et ail (]2008T ) 
should be bright in the FUV. We obtained FUV obser- 
vations using the solar blind channel (SBC) of the Ad- 
vanced Camera for Surveys (ACS) on the Hubble Space 
Telescope (HST) to search for the predicted FUV emis- 
sion. We describe the observations and analysis in §2. 
We discuss the implications in §3. 

2. OBSERVATIONS 

Observations were performed on 2008 May 01 using the 
ACS/SBC on HST centered on the position of the X-ray 
source in MF 16 (GO program 11200, PI: Kaaret). The 
spectral element used was the F140LP long-pass filter 
which covers the 1350-2000A band. This filter was cho- 
sen because it has a much smaller background sky rate, 
blocking out geocoronal Lyman a emission, than the fil- 
ters that extend to shorter wavelengths. The exposure 
time was 2760 seconds. 

We also extracted several images of the field from the 
HST archives, these include an ACS Wide Field Camera 
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Fig. 1.— Images of MF 16. The panels show A: FUV (F140LP filter), B: B-band (F450W filter), C: V-band (F555W filter), and D: 
I-band (F814W filter). The small circles mark the counterparts to the X-ray source and have radii of 0.15". The ellipses mark the extent 
of the nebula with a major diameter of 2.0" and a minor diameter of 1.4". The arrow in panel A points North and is 1" long. 



(WFC) image in the F814W (I-band) filters obtained on 
2004 July 29 and the Wide Field and Planetary Camera 
2 (WFPC2) images in the F450W (B-band), F555W (V- 
band), and F814W (I-band) filters all obtained on 2001 
June 8. For the ACS data, we used the dithered images 
directly from the archive. For the WFPC2 data, we com- 
bined the 4 images obtained in each filter and performed 
cosmic ray rejection using multidrizzle in Pyraf. 

We re-analyzed an archival observation made with the 
Advanced CCD Imaging Spectrometer (ACIS) on the 
Chandra X-ray observatory on 2001 September 07 UT 
(ObsID 1043). This observation was chosen because it 
had the longest exposure time, 59 ks, of those avail- 
able in the archive, allowing the deepest search for X- 
ray counterparts to optical sources as needed for the X- 
ray/optical image alignment. The Chandra data were 
subjected to standard calibration using CIAO v4.0: ap- 
plying the aspect solution, generating an instrument map 
weighted by an absorbed power law, and generating an 
exposure map. There were no background flares, so the 
entire exposure time was used. 

2.1. Astrometry 

The Two-Micron All Sky Survey 

f2MASS ISkrutskie et all l200l catalog was chosen 
for absolute astrometry correction over other catalogs, 
e.g., USNO-B1.0, because of its superior accuracy. 
Because the SBC has a small (34"6 x 30'.'8) field of view, 
there are no 2MASS stars on our image. Instead, abso- 
lute astrometry was performed using the WFC/F814W 
image, and then the SBC and Chandra images were 
aligned to the WFC/F814W image. 

Using the Graphical Astronomy and Image Analysis 
(GAIA) tool, we detected stars on the WFC/F814W im- 
age and then searched for point-like counterparts with 



TABLE 1 
Photometry 



Filter 


Count rate 


Flux for Ay 


= 1 


54 




(c/s) 


(erg cm -2 s — 


l A- 


- 1 ) 


F140LP (FUV) 


0.24 ±0.04 


(4.8 ± 0.7) x 


10" 


16 


F450W (B) 


0.41 ±0.06 


(2.8 ± 0.4) x 


io- 


17 


F555W (V) 


0.69 ±0.10 


(1.3 ± 0.2) x 


io- 


17 


F814W (I) 


0.33 ±0.05 


(2.7 ±0.4) x 


10" 


18 



magnitudes K < 15 in the 2MASS catalog, finding 20 
matches. These matches were used to compute a new 
aspect solution for the F814W. The fit rms was 1.9 pix- 
els or 0.9". Eighty- nine point sources were found in the 
SBC image, of which 11 are also appear in the F814W 
image. These 11 matches were used to transfer the aspect 
solution to the SBC image. 

We found 40 X-ray sources within the field of the 
F814W image by running the wavdetect tool in CIAO 
on a 0.3-8 keV band image. Of these 40, 3 have coun- 
terparts in the WFC image (excluding the X-ray source 
in MF 16). However, one of the objects is extended and 
another is faint and has a large point-spread function due 
to its off-axis position in the Chandra image. So the one 
bright, point source was used to correct the astrometry. 
First, the positions found from wavdetect were corrected 
by subtracting the systematic error found by running a 
Marx simulation of the Chandra telescope. The Chan- 
dra image was then translated, by 0.2" in right ascension 
and —0.3" in declination, to match the F814W image. 
The corrected position of the X-ray source in MF 16 is 
R.A.=20:35:00.752, decl.=+60:ll:30.55 (J2000). From 
the wavdetect position errors, we estimate the 90% er- 
ror radius to be 0.4". 



ULUV source in NGC 6946 



3 



2.2. Photometry 

Figure Q] shows images of MF 16 in the FUV (F140LP 
filter), B-band (F450W filter), V-band (F555W filter), 
and I-band (F814W filter). The nebula and a counter- 
part to the X-ray source are clearly visible in each band. 
The nebula appears extended along the East- West direc- 
tion. The ellipse drawn on the figure has a major diame- 
ter of 2.0" and a minor diameter of 1.4" corresponding to 
a physical size of 50 p c by 35 pc at a distance of 5.1 Mpc 
(|deVaucouleurslll979l) . 

We performed aperture photometry on the X-ray 
source counterparts, extracting counts from a circle with 
a radius of 0.15". We use the WFPC2 F814W image, 
rather than the ACS/WFC one, since the source may 
be variable and the three WFPC2 images were obtained 
over a span of less than 5 hours. We subtracted the 
nebular flux, scaled by area, from the point source flux. 
If, instead, we use a background region external to the 
nebula, i.e. retain the nebular contribution within the 
0.15" extraction circle, the FUV counts increase by 11% 
and the BVI counts increase by ~ 30%. We relied on 
the synphot package in IRAF for the aperture correc- 
tion for the ACS data. For the WFPC2, we determined 
the aperature correction (from a radius of 0.15" to 1.0") 
using the star USNO B1.0 1501-0283184 also located on 
the planetary camera chip. The same star was also used 
to check our photometry in the B and I bands. Nebular- 
background subtracted count rates for the point source 
for each continuum filter are given in Table [TJ We esti- 
mate a 15% error on the source count rates, primarily due 
to uncertainties in the background subtraction (due to 
the presence of the nebula). The background-subtracted 
count rate in the F140LP filter for the nebula, excluding 
the point source, is 0.78 c/s. For the point source, we 
find equivalent magnitudes of B = 22.73 and V = 22.70. 
The V magni tude is in good agreement with V — 22. 64 
as quoted by I Blair. Fesen. fc Schlegel (|2001| ). while our 
B-magnitude is somewhat brighter than their value of 
B = 23.10. 

We modeled the intrinsic source flux as a power-law, 
F\ oc A -2 and used the calcphot command in the syn- 
phot package to relate the measured count rates to in- 
trinsic source fluxes. We note that the inferred fluxes 
change by less than 1% if we change the power-law spec- 
tral index to F\ oc A -4 . The filter pivot wavelengths 
are 1533 A for F140LP, 4557 A for F450W, 5443 A 
for F555W, and 7996 A for F814W. The Galactic in- 
terstellar extinctio n along the line of sight to MF16 is 
A v = 1.14 mag (jSchlegel et all 11998ft . Reported val- 
ues of the additional, intrinsic extinction, mostly found 
by measuring decrements between Balmer lines in the 
nebular spectrum, vary from 0.2 to 0.66 magnitudes 
(iBlair. Fesen. fc Schlegdll200H IRo berts fc Colbertl 120031: 
lAbolmasov et al.l I2008D . Following lAbolmasov etlti! 
( 20081) . we adopt an exti nction of Ay = 1.54 mag 
and th e extinction curve of ICardelli. Clayton, fc Mathisl 
(1989). The point source fluxes are given in Table [1] 
for Ay — 1.54 mag. The l ower extinction value dis- 
cussed by lAbolmasov et al.l (|2008l ). Ay = 1.34 mag, 
decreases the intrinsic flux by 41% in the FUV, 22% 
in B, 18% in V, and 11% in I. The lower extinction 
may be valid for the point source, if source preferen- 
tially ionizes the nearby regions of the nebula. Adopt- 
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Fig. 2. — Spectral energy distribution of the point-like source 
at the center of MF 16. The dereddcncd fluxes from Table [T] are 
plotted as diamonds. The dashed line is a power-law with F\ oc 
A~ 41 . The dotted curve is the sum of two blackbody spectra with 
temperatures of 140,000 K and 26,000 K. The solid curves show 
multicolor disk blackbody spectra. The curve with higher flux at 
200 A is for p = 0.75 and the lower one is for p = 0.5. 

ing an extinction at the high end of the reported val- 
ues, A v = 1.8 mag jBlair. Fesen. fc Schlegell 120011: 
IRoberts fc Colberdl2003T) . increases the intrinsic flux by 
84% in the FUV, 33% in B, 26% in V, and 15% in I. 

3. DISCUSSION 

FUV imaging of MF 16 reveals both a point-like source 
coincident with the X-ray source and a surrounding neb- 
ula with an angular extent of 2.0" by 1.4". The point 
source has a flux of 5 x 10 -16 ergs -1 cm -2 A -1 and the 
nebula has a flux of 1.6 x 10 -15 ergs -1 cm -2 A -1 , both 
quoted at 1533 A, corresponding to the pivot wave- 
length of the F140LP filter and assuming an extinction 
of Ay = 1.54. The monochromatic luminosity estimate 
is XL\ = 1.0 x 10 40 ergs -1 , suggesting the presence of an 
ultraluminous UV source ( ULUV) . The mea s ured fl ux is 
close to that predicted by lAbolmasov et all (|2008fi near 
1000 A. 

Figure [2] shows the spectral energy distribution of the 
point-like source at the center of MF 16 using the data 
from the FUV (F140LP), B (F450W), V (F555W), and I 
(F814W) bands. The dereddened fluxes are from Table[TJ 
We note that the BVI data are quasi-simultaneous, ob- 
tained within 5 hours, while the FUV data were obtained 
7 years later. Thus, the non-simultaneity of the observa- 
tions should be o f conc ern in fitting the data. However, 
iFridriksson et all (|2008f ) report that the rms variation in 
X-ray flux of this source in Chandra observations spread 
over several years is only 8%. Also, the I-band flux from 
the ACS/WFC measured in 2004 agrees with that mea- 
sured with the WFPC2 in 2001 within 25%. Thus, fitting 
of these non-simultaneous data appears justified. 

A fit of a power-law, F x = A(A/6000A) Q to the BVI 
data dereddened assuming Ay = 1.54 gives best fit 
parameters a = -4.1 ± 0.4 and A = (8.8 ± 0.8) x 
1Q- 18 ergs" 1 cm" 2 A" 1 , see Fig. gj Using Ay = 1.8 
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makes the curve somewhat steeper, a — —4.4 ± 0.4, 
while Ay — 1.34 leads to a slightly shallower slope, 
a = —3.9 ± 0.4. In all good fit is obtained with 

X 2 < 1. The exponent of the the power- law is consis- 
tent with that expected for the Rayleigh- Jeans tail of a 
blackbody spectrum, a = —4. 

We fit all four data points with a single blackbody 
spectrum. This produces a good fit with a temper- 
ature T = 31000 K and a bolometric luminosity of 
3.4 x 10 39 ergs -1 . However, due to the relatively low 
temperature, the spectrum produces an insignificant flux 
of radiation capable of ionizing He. A second source pro- 
ducing radiation at shorter wavelengths is required to 
produce the Hen A4686 luminosity seen from the nebula. 
Thus, we added a hotter blackbody with T =140,000 K 
and L = 1.2 x 10 4 ergs~ 1 as suggested b y the emis- 
sion line modeling (|Abolmasov et all 120081 ). The fit is 
acceptable, x 2 /DoF = 3.8/2, and the cooler blackbody 
temperature is reduced to T =26000 K and its luminos- 
ity to 2.2 x 10 39 ergs -1 . These values are for Ay = 1.54. 
Varying the extinction in the range Ay = 1.34 — 1.8 and 
the hot blackbody parameters w ithin the certainty range 
given bv lAbolmasov et al.l (|2008t ). allows temperatures in 
the range 23,000-32,000 K and luminosities in the range 
1.1 — 4.6 x 10 39 ergs -1 . These properties are consistent 
with a high mass star, specifically a late O or early B 
(07 to B3) supergiant, that could be the companion star 
in the binary system. It would be of interest to obtain 
an optical and/or UV spectrum of the point source to 
search for spectral lines that would permit classification 
of the star, or rule out interpretation in terms of a stellar 
companion. 

The spectra of accretion disks around black holes are 
often described by a multicolor disk blackbody (MCD) 
spectrum where the local effective temperature is a de- 
creasing function of radius, usually T(r) oc r~ p , be- 
tween an inner disk radius, Ri n , where the tempera- 
ture is T in and an outer radius, i? ou t, with tempera- 
ture T out . At long wavelengths, hc/X <C fcT out , the 
spectrum has the form of a Rayleigh- Jeans tail and 
at intermediate wavelengths, fcT out <C hc/X <C fcTi n , 
the spectrum flattens. The observed FUV/optical spec- 
trum of the point source is compatible with this form 
if the outer disk temperature lies between the FUV and 
B bands. The FUV/optical data can constrain MCD 
model spectra only in regards to the outer disk temper- 
ature and the flux in the FUV/optical range. A con- 



straint on the short wavelength end of the spectrum 
can be obtained from the ionizing luminosity needed to 
produce the nebular Hell A4686 emission. The solid 
lines in Fig. [2] show MCD spectra consistent with the 
FUV/optical data and with a UV luminosity of i(20 < 
A < 228A) w 4 x 10 39 ergs 1 as needed t o produce the 
observed Hen line (Abolmas ov et al1l2008f) . We consider 
two disk temperature profiles: a standard profile with 
p = 0.75 and a profile with p = 0.5 as found in spec - 
tral fits to a number of ULXs (|Vierdavanti et al~l f2006y 
Both fits are acceptable with x 2 < 3.4. For p — 0.75, 
the bolometric luminosity is LboI = 1-2 x 10 ergs -1 
and T out = 16000 ± 2000 K. However, the disk is un- 
usually narrow, -R ut/-Rin ~ 40. For p = 0.5, we find 
that X B oi = 8 x 10 40 ergs -1 , T out = 21000 ± 2000 K, and 
-Rout /Rin ~ 6000. The true situation may be more com- 
plex with p ~ 0.5 only in the inner disk and R ou t/Rm 
intermediate between these two extremes. It is also be 
possible to fit the data using models where the outer disk 
is irradiated. However, the complexity of irradiation, and 
the number of parameters in such models, preclude ob- 
taining useful constraints with the available data. In all 
cases, the outer disk radius is of order 10 12 cm. 

In summary, two possible scenarios to explain the ob- 
servational data are that the emission is the sum of a 
hot blackbody from the compact object and stellar emis- 
sion from the companion star or that the emission is due 
to a multicolor disk with an outer disk radius ~ 10 12 cm. 
The small number of broad band measurements currently 
available is inadequate to distinguish between these two 
cases. Additional observations are needed to determine 
the physical nature of the source. 
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